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ABSTRACT 

The  density  and  parallel  temperature  of  an  electron  beam 
passing  through  a  magnetic  guide  field  and  a  linear  wiggler  are 
determined  as  functions  of  distance  from  the  anode.  Starting  with 
a  given  beam  emittance  and  phase  space  distribution  at  the  anode, 
single-particle  orbits  are  employed  to  propagate  the  beam  along 
the  drift  tube.  Finite  gyro-orbit  effects  induce  substantial  modulati 
in  density.  Additionally,  the  evolution  of  an  initial  spread  in 
the  parallel  speed  of  the  electrons  is  sensitively  dependent  on 
the  wiggler  strength  and  contributes  to  further  density  fluctuation. 
The  beam's  parallel  temperature,  a  critical  factor  in  determining 
laser  efficiency  and  gain,  is  found  to  increase  significantly  when 
the  wiggler  strength  is  raised. 
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I.  INTRODUCTION 


In  a  free-electron  laser,  a  beam  of  electrons  of  relativistic  mass 
factor  y  is  passed  through  a  magnetic  wiggler  field  of  wavelength  X  resulting 
in  emission  of  electromagnetic  waves  of  wavelength  Xracj  s  X/(2y  ]  whose 
intensity  is  dependent  on  the  beam  density.  For  the  most  part,  previous 
theoretical  analyses  of  the  free-electron  laser  instability  have  employed 

% 

highly  simplified  models  which  assume  this  density  is  spatially  uniform. 
Noting  that  any  function  of  the  constants  of  motion  of  a  single  electron 
is  a  steady  state  solution  of  the  Vlasov  equation,  Davidson  and  Uhm^  were 
the  first  to  construct  self-consistent  equilibria  for  electron  beams  propa¬ 
gated  in  an  axial  guide  field  and  a  helical  wiggler.  The  resulting  spatial 
modulation  of  the  electron  density  was  determined  for  different  choices 
of  the  functional  dependence  of  the  distribution  function  on  the  constants 
of  motion. 

In  this  paper,  we  examine  the  experimental  situation  in  which  the 
directional  flux  of  electrons  (and  beam  emittance)  at  the  anode  is  specified 
and  the  beam  is  allowed  to  propagate  through  a  drift  tube  immersed  in 
a  strong  axial  guide  field  and  a  linear  (i.e.,  planar)  wiggler.  Using 
nonlinear  particle  orbits  obtained  earlier,  the  distribution  function 
of  the  electrons  at  any  point  along  the  drift  tube  is  determined  from 
the  equation  of  continuity  in  phase  space.  The  electron  density  is  then 
obtained  by  an  integration  over  momenta.  Unlike  the  calculations  in  Ref. 

1,  an  analysis  using  these  orbits  permits  us  to  treat  the  case  where  there 
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is  a  resonance  between  the  gyromotion  of  the  electrons  and 
axial  periodicity  in  parallel  particle  velocity  induced  by 
In  addition,  the  radial  dependence  of  the  wiggler  field  is 


the  effective 
the  wiggler. 
allowed  for. 
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II.  MATHEMATICAL  ANALYSIS 


We  begin  by  choosing  a  rectangular  coordinate  system  with  the  anode 
at  its  origin  and  the  drift  tube  oriented  along  the  positive  z  axis.  In 
Ref.  2  it  is  argued  that  in  the  presence  of  an  axial  guide  field  it  is 
natural  to  introduce  the  variables  y,  0,  X,  Y  which  are  defined  as  follows 


y  =  (Zmfto)"1  [(Py  +  m  tt0y)2  +  p2]  ' 

A  y 

e  =  arctan  [(Px  +  m  ftQy)/py] 

_  i 

X  =  x  +  (m  fi0)  py 
Y  =  -(mn0)'1Px  . 


(1) 


Here  f20  =  qH0/mc  is  the  gyrofrequency  of  a  particle  with  charge 

q  and  (rest)  mass  m  in  a  magnetic  field  of  intensity  H0.  The  canonical 

momenta  P  ,  p  ,  p  are  related  to  the  particle's  velocity  components 
x  y  z 

vy,  v  ,  v  through  the  expressions 
X  y  z 


px  =vmx  +  (q/c)Ax  =  px  +  (p/c)Ax 


py,z  =^y>z 


A  =  (Ax,Ay,Az)  =  [-H0y  +  6Hk 


sin(kz)cosh(ky) ,0,0] 


where 
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is  the  vector  potential  of  a  linear  magnetic  wiggler  of  the  form 
[0,  6Hcos(kz)cosh(ky) ,  -  SHsin(kz)sinh(ky)]  with  wavelength  2Trk-1  along 
the  z-axis,  superimposed  on  a  guide  field  (0,0,Ho).  Note  that  in  the 
absence  of  the  wiggler  field  (SH  -0)  the  foregoing  definitions  of 
y,0,X,Y  reduce  to  the  usual  expressions  for  the  magnetic  moment,  gyroangle, 
and  the  x  and  y  coordinates  of  the  guiding  center. 

It  is  shown  in  Ref.  2  that  for  sufficiently  small  wiggler  fields 
(measured  by  e  =  6H/H0)  the  motion  of  the  electrons  can  be  locked  onto 
a  single  resonance,  with  the  particles  strongly  influenced  by  this  resonance 
but  mildly  so  by  adjacent  harmonics.  This  resonance  condition  has  the 
form  &e+kz s  0  for  some  integer  £,  where  z  is  the  particle's  axial  velocity. 
It  is  further  shown  that  in  this  case  the  motion  is  nonstochastic,  and 
explicit  integration  of  the  equations  of  motion  is  possible.  For  free- 
electron  laser  applications  it  is  imperative  to  limit  the  strength  of 
the  wiggler  in  order  to  avoid  stochastic  motion. 

O  A 

In  the  experiment  conducted  by  Roberson  et  al.  *  at  the  Naval  Research 
Laboratory,  which  provided  the  original  motivation  for  this  work,  a  current 
Ig  of  electrons  (=500  A)  is  accelerated  through  a  potential  difference 
of  350  kV  (y  =  1.68)  before  impinging  on  an  anode  of  radius  R0  =  0.5  cm 
and  entering  a  60  cm  drift  tube.  From  these  numbers  it  follows  that 

'i  -  f~- 2  * 10-2 

0  0 

where  1^  is  the  Alfven  current,  v  =  N(e  /me  )  is  Budker's  parameter, 

N  =  2tt  /  *°  rn(r)dr  is  the  number  of  electrons  per  unit  axial  length,  and 
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n(x,y,z)  is  the  electron  density.  With  v/y  «  1,  electric  and  magnetic 

5 

self-field  effects  are  measured  by  the  dimensionless  parameter 
n  =  (2y<jo  2/ft2)[l-(v  /c)2]  where 

Ho  ^ 

o)p  =  (4Trne2/m)5  (2) 

is  the  electron  plasma  frequency.  For  the  NRL  experiment  q  =  1/2.  Further, 
in  the  experiment  the  magnetic  field  is  allowed  to  build  to  its  full  value 
in  a  transition  region  preceding  the  drift  tube  entrance.  Both  the  work 
described  here  and  in  Ref.  2  is  idealized  in  that  neither  self-field  effects 
nor  the  presence  of  a  transition  region  are  allowed  for.  In  addition, 
because  of  uncertainties  in  the  magnitude  of  various  experimental  parameters 
we  shall,  in  what  follows,  estimate  appropriate  values  when  necessary. 

The  electron  density  at  r  is  given  by 

n(r)  =  /d3p  f(r,p)  (3) 

where  f(r,p),  the  electron  distribution  function,  satisfies  the  conservation 
equation  in  six-dimensional  phase  space.  Given  this,  it  is  clear  that 
the  electron  density  as  a  function  of  the  distance  along  the  drift  tube 
in  the  steady  state  may  be  written 

n(?)  =  /d3  p  ?(?,?)]  (4) 

where  [^(r,p),  ?(r,p)]  is  the  phase  space  orbi t  of  an  electron  which  starts 
at  the  anode  and  reaches  the  point  (r,p)  =  (x,y,z,p)  at  the  axial  location 
z. 
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The  final  assumption  of  this  work  involves  specifying  the  distribution 
function  f  at  the  anode.  We  shall  assume 

f  (anode)  =  C0(Ro-R)6[H-(P^+Py+p^)]exp(-<)>)  (5) 

where  C  is  a  normalization  constant,  ©(•)  is  the  Heaviside  unit  step  function, 

R0  is  the  anode  radius,  R  e  (x  +y  )2  is  the  radial  distance  of  the  point 

(x,y,z)  from  the  axis  of  the  drift  tube,  H  is  related  to  the  Hamiltonian 
2 

H  =  ymc  by 

it  -  H*,C2-m2C2 

=  [Px+mft0y-m6fik-1sin(kz)cosh(ky)]^  +  p^+p^ 
with  SQ  E.q5H/mc,  and  cf>  is  given  by 


4> 


+ 


(6) 


In  (6),  pl(  and  p±  are  the  spread  in  the  parallel  and  perpendicular 
momenta  of  the  electrons  emerging  from  the  anode,  with  the  parallel  momenta 
centered  about  the  resonant  value  pfa  =  £m|ft0|/k  (Ref.  2).  The  Heaviside 
function  in  the  expression  for  f  indicates  that  the  beam  entering  the 
drift  tube  has  radius  R0  and  a  uniform  density  over  the  anode  area.  The 
6-function  indicates  that  the  beam  is  monoenergetic,  with  a  fixed  y. 

The  experimentally  determined  spread  in  the  perpendicular  speed  is  roughly 
given  by  (p^/mc)  =  0.08  -  0.24;  by  a  simple  kinematical  argument  one  finds 
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2 

that  p(|  =  1/2 (p ^/Pj^)  s  whence  the  spread  in  parallel  speed  is 
determined. 

Given  the  distribution  function  and  the  (canonical)  transformation 
(1),  we  proceed  with  the  determination  of  n(r).  Since  the  Jacobian  of 
(1)  is  m|fi0|  (so  that  dpxdpydxdy  =  m|fi0 |dOdydXdY)  we  insert  f  into 
expression  (4),  integrate  over  pz,  and  obtain 

m|ft0|  i 

n(r)  =  C~2 -  /dedydXdY0{Ro-R[£(r,p),  ?(r, p)]}  p~X 

(7) 

•  exp{-<}>[^(r  ,p) ,  ?(r,p)]}  <5  (x-X-pcos0)6(y-Y+psi  n0 ) , 

where  the  factor  1/2  arises  from  the  neglect  of  backward  propagating 
particles,  and  the  gyroradius  p  =  (2y/m£20r  relates  the  particle  (x,y)  and 
guiding-center  (X,Y)  coordinates  via 

X  =  x  -  pcose 
Y  =  y  +  psine  . 


In  order  to  complete  the  integration  in  (7)  an  explicit  description 
of  the  particle  orbit  [i<(r, pj,  £(r,p)]  is  needed.  The  calculation  of 
these  orbits  is  presented  in  the  Appendix.  Having  determined  them,  (7) 
reduces  to  a  double  integral  in  y  and  0  which  can  be  evaluated  using 
numerical  techniques.  Although  y  appears  to  vary  over  a  semi-infinite 
interval,  the  presence  of  a  Gaussian-like  integrand  effectively  reduces 
the  region  of  integration  to  finite  size.  As  z  is  varied,  the  size  of 
this  region  is  adjusted  dynamically  in  order  to  include  all  significant 
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contributions  of  the  integrand.  Various  mesh  densities  were  employed 
in  order  to  ensure  convergence. 

Having  computed  n(r)  it  is  straightforward  to  determine  the  parallel 
temperature  Tz  given  by 

Tz(r)  e  I  /d3pf(?,f)-fin[^r  -  V7(?)]2 

^  ii  Ym  z 


where  the  mean  parallel  velocity  Vz  is 


The  quantity  Tz  is  a  measure  of  the  spread  in  parallel  velocities  in  the 
beam  frame.  In  order  that  the  electromagnetic  radiation  emitted  by  the 
electrons  be  as  coherent  as  possible,  and  to  allow  operation  in  the  col¬ 
lective  (i.e. ,  Raman)  regime,  Tz  must  be  small. 

Profiles  of  the  density  and  parallel  temperature  for  various  wiggler 
field  strengths  and  harmonics  (given  by  e  and  i)  are  presented  in  the 
figures.  These  are  reviewed  in  the  next  section.  Plotted  results  are 
accurate  to  1%. 
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III.  DISCUSSION 


In  Fig.  1,  the  beam  density  (in  arbitrary  units)  is  plotted  as  a 
function  of  the  axial  distance  z  from  the  anode  along  the  drift  tube  for 
e  =  1/1024.  For  the  upper  curve  the  observation  point  at  which  the  density 
is  measured  has  coordinates  x  =  0,  y  =  0.5  while  the  lower  curve  has 
x  =  0,  y  =  1.0  x  10  .  All  linear  dimensions  are  given  in  centimeters. 

The  substantial  density  modulation  evident  in  the  upper  curve  can  be  under¬ 
stood  by  examining  the  nature  of  the  particle  orbits  using  Fig.  2.  In 
this  diagram,  circle  A  represents  the  anode  (chosen  to  have  radius  1.2) 
and  0  is  the  observation  point  (coordinates  x  =  0,  y  =  0.5).  (We  use 
a  somewhat  larger  anode  radius  in  our  numerical  work  in  order  to  compensate 
for  uncertainties  in  the  experimental  values  of  y  and  the  spread  in  v,). 
Consider  gyro-orbits  which  enclose  circle  R,  pass  through  observation 
point  0  and  have  their  guiding  center  on  the  y-axis  (circle  R  for 

X/ 

example).  Particles  in  such  orbits  necessarily  move  beyond  the  anode 
radius  during  part  of  their  trajectory.  Since  the  beam  has  little  spread 
in  parallel  momentum  by  far  the  greatest  number  of  particles  have  gyroradii 
equal  to  0.453,  that  of  particles  exactly  at  resonance  ( p  =p.  =£m|£20 1 /k) . 
Circles  R  and  Ru  are  gyro-oribts  that  bound  this  most  probable  value, 
and  are  populated  with  ^10”^  of  the  number  of  particles  at  resonance. 

As  we  move  forward  in  space  from  the  anode,  the  particles  in  the  relevant 
orbits  (denoted  by  the  shaded  regions)  move  outside  A  and  the  density 
at  0  falls.  For£=  2,  the  resonance  condition  JlS+kzsO  implies  that  in 
one  wiggler  wavelength  (3  cm)  the  particles  complete  about  half  of  their 
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gyro-orbit  and  the  density  should  be  minimum.  Over  a  distance  of  two 
wiggler  wavelengths  the  particles  return  to  their  original  locations  (in 
the  limit  where  all  parallel  momenta  are  identical)  and  the  density  is 
again  at  its  peak. 

Because  a  small  spread  in  parallel  momenta  actually  exists,  particles 
completing  their  orbits  return  to  their  original  positions  at  slightly 
different  times,  producing  a  smaller  modulation  of  density  evident  in 
the  shorter  wavelength  oscillations  in  the  upper  curve  of  Fig.  1  and  the 
latter  part  of  the  lower  curve. 

Of  course,  there  are  other  regions  besides  the  shaded  area  in  Fig. 

2  which  contribute  to  the  plotted  density  modulations.  We  have  simplified 
our  discussion  and  diagram  to  keep  the  argument  clear.  By  recalling, 
however,  that  the  only  significant  gyroradii  are  approximately  equal  to 
0.453  and  that  the  observation  point  is  0,  it  is  easy  to  see  that  the 
only  particles  contributing  to  the  density  modulation  lie  in  an  appropriately 
expanded  shaded  region  symmetrically  placed  about  the  y-axis.  For  instance, 
since  the  radius  of  circle  L  in  Fig.  2  is  0.85,  practially  no  particle 
with  a  guiding  center  lying  on  the  y-axis  to  the  left  of  0  whose  orbit 
passes  through  0  moves  out  of  A.  Such  a  particle,  therefore,  does  not 
contribute  to  the  observed  density  depletion. 

The  wiggler  acts  to  alter  the  speed  of  the  particles  as  they  move 
along  the  drift  tube.  By  increasing  the  amplitude  of  the  wiggler  field 
an  enhancement  of  the  shorter-wavelength  modulation  seen  in  Fig.  1  is 
expected.  This  is  clearly  evident  in  Figs.  3  and  4  in  which  e  =  1/128 
and  e  =  1/64  respectively.  It  is  this  short-wavelength  modulation  that 
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is  apparently  discussed  in  Ref.  1,  for  in  that  work  it  is  assumed  there 
is  no  spread  in  the  parallel  momenta  of  the  electrons.  We  further  note 
that,  since  Xracj  is  typically  much  smaller  than  the  wiggler  wavelength, 
the  observed  modulations  in  density  not  only  should  modify  the  gain  but 
may  also  adversely  affect  the  coherence  and  propagation  characteristics 
of  the  emitted  radiation. 

2  o 

In  Fig.  5  the  normalized  parallel  temperature  T  =  (ck/£fi0)  (yT  /me  ) 
is  plotted  as  a  function  of  axial  distance  z  for  e  =  1/1024  and  e  =  1/128.  In 
both  cases  the  observation  point  has  coordinates  (x,y)  =  (0,0.5).  One 
discerns  the  same  six-centimeter  modulation  evident  in  the  density  plots. 
Again,  these  fluctuations  are  due  to  the  spread  in  the  parallel  speeds  of 
the  particles  and  the  fact  that  some  particles  periodically  vacate  the  ob¬ 
servation  point  and  move  outside  of  the  nominal  beam  radius.  The 
enhancement  of  the  spread  in  parallel  speeds  caused  by  the  wiggler  is 
apparent  by  comparing  the  two  curves.  Note  that,  although  the  upper 
curve  corresponds  to  a  wiggler  almost  an  order  of  magnitude  stronger  than 
the  lower  curve,  both  wiggler  fields  are,  nonetheless,  relatively  weak. 
Furthermore,  while  the  parallel  temperature  is  seen  to  increase  by  about 
two  orders  of  magnitude  it  is  still  not  very  large  in  absolute  terms, 

roughly  10  keV  at  z  =  24  cm,  for  e  =  1/128. 

6  7 

As  is  well-known,  *  the  high  gain  and  conversion  efficiency  associ¬ 
ated  with  the  operation  of  a  free-el ectron  laser  in  the  Raman  regime  is 
predicated  upon  the  availability  of  a  cold  electron  beam  having 
( 6yz ) /yz  <  (A/2y3/2)(cop/2c) ,  which  may  be  rewritten  as 
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(Svz) 


0)„  ^  2 
fa2  )(#-) 


v  ^  3/2  vco  /  'v 

z  y  '  rad  z 


(8) 


Here  wracj  is  the  frequency  of  the  emitted  radiation,  Wp  is  defined  in  (2),  and 

2  2  -x 

6yz  and  6vz  denote  the  spread  in  yz  =  (l-Vz/c  )  2  and  vz,  respectively. 

2  2 

With  (<Svz)/yz)  =  Tz/ ( ymvz ) ,  expression  (8)  can  be  rewritten  in  terms  of 

4  2-1 

the  normalized  parallel  temperature  as  T  =  (c/vz)  Y  •  ^or 

3  4 

the  NRL  parameters  ’  T  s  0.064.  Thus,  for  the  two  cases  depicted  in 
Fig.  5  the  beam  is  sufficiently  cold  to  allow  operation  in  the  Raman 
regime. 

Figure  6  shows  how  the  density  modulations  are  affected  by  an  increase 
in  the  parallel  momentum  of  the  particles.  With  the  l  =  3  harmonic  resonant, 
y  =  2.17  and  e  =  1/128  a  typical  gyroradius  is  the  same  as  in  the  £  =  2 
case.  For  the  upper  curve  in  Fig.  6  the  depth  of  modulation,  defined 

as  (nmax  -  nroJ/("max  +  nmin>’  1s  ver^  c1ose  t0  that  of  F19-  3-  ™s 
is  to  be  expected  since  the  gyroradius  and  anode  radius  are  the  same  in 

both  cases.  Since  £  =  3,  one  also  expects  the  density  modulations  to 

have  a  9  cm  wavelength,  which  is  clearly  evident  in  Fig.  6. 
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IV.  CONCLUSION 


The  successful  operation  of  a  free-electron  laser  is  strongly  de¬ 
pendent  on  the  quality  of  the  electron  beam  employed.  It  is  known  that 
both  the  laser  gain  and  the  coherence  of  the  radiated  electromagnetic 
energy  may  be  adversely  affected  by  density  fluctuations  in  the  beam  or 
spreads  in  parallel  electron  speed.  In  this  paper  we  have  examined  an 
electron  beam  propagated  through  an  axial  guide  field  and  a  linear  (i.e., 
planar)  wiggler.  By  employing  single  particle  orbits,  the  density  and 
temperature  of  the  beam  as  functions  of  distance  from  the  anode  have  been 
determined.  We  present  results  which  indicate  that,  unless  care  is  taken, 
rather  substantial  modulations  in  beam  density  can  be  expected.  In  examining 
the  density  both  on  and  off  axis  for  various  wiggler  field  strengths, 
we  have  found  that  these  modulations  increase  both  as  the  edge  of  the 
beam  is  approached,  and  as  the  wiggler  strength  is  increased.  This  occurs 
for  wigglers  whose  field  strength  is  relatively  small  compared  to  the 
guide  field.  As  a  rough  criterion  for  suppressing  significant  density 
fluctuations,  we  conclude  that  the  electron  gyroradius  must  be  a  small 
fraction  of  the  beam  radius.  Moreover,  the  initial  spread  in  the  parallel 
speeds  of  the  electrons  must  be  small  and  the  variation  in  this  parallel 
speed  induced  by  the  wiggler  must  be  comparable  to  this  spread  or  smaller. 
This  may  result  in  the  need  to  limit  the  strength  of  the  wiggler,  otherwise 
operation  in  the  Raman  regime  may  not  be  possible. 
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APPENDIX:  SINGLE-PARTICLE  ORBITS 


In  this  appendix,  the  single-particle  equations  of  motion  are  inte¬ 
grated  and  the  particle  orbit  expressed  as  a  function  of  the  distance 
from  the  anode. 

It  is  shown  in  Ref.  2  that  if  the  wiggler  field  strength  is  small 
enough,  so  that  the  particle  motion  is  influenced  primarily  by  the  5,-th 
resonance,  then  the  equations  of  motion  reduce  to 


(Al) 


where 


2 

=  -  (I+J)  , 


k 


(A2) 


and 


Here,  the  effective  potential  V  ^  is  given  by 


o  2  o  d  I  n  p 

WI>=<H*-I  >  -  [7 MkY>  f -a r] 
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with 


C  =  kp  =  £[2(I+J)] 


1 


(A3) 


F£(x)  =  i[exp(x)  -  (-if  exp(-x)]  , 

and  1^  =  I.  (5)  the  modified  Bessel  function  of  the  first  kind  of  order  i  . 

Since  the  resonant  value  of  the  axial  momentum  equals  -£mf20/k  (note 

fto  <  0  for  electrons)  it  follows  from  (A2)  that  1=0  at  exact  resonance. 

In  Ref.  2  several  possible  forms  for  the  effective  potential  are 
sketched.  For  bounded  motion  of  amplitude  a  and  frequency  a  in  a 


potential  well  having  a  minimum  at  I  ,  (Al)  may  be  approximately  integrated 
to  obtain 


I  =  I  +  a  si n[a(t-x)+n]  , 


(A4) 


where 


(A5) 


•• 


sin(n)  =  ^ 


with  the  boundary  condition  that  I  =  I  at  t  =  x.  (In  our  numerical 
work  bounded  potential  wells  are  approximated  by  parabolas.  This  permits 
explicit  determination  of  expressions  (A4)  and  (A5)).  Inserting  (A4) 
into  (A2)  and  using 
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p  =  ym 


dz 

dx 


we  obtain 


z  =  z  (1-1  )(t-x)  -  ^g-a{cos[a(t-x)+n]-cos(n)  }  . 


SUQ 


(A6) 


The  equation  of  motion  for  the  x- component  of  the  guiding  center  may 
2 

be  expressed  as 


X  = 


-Ao  F£+l^kY^  /LI  t2\ 

2^  ~TW  (H£  '  1  ] 


which,  upon  inserting  (A4)  and  integrating,  yields 


„2  F  -i  (kY)  2  9a: 

X  =  X-^-T^kYT<^-  I-V)(t-x)+^{cos[a(t-T)+n] 


cos(n)} 


+  ■§£  (sin[2a(t-x)+2n]  -  s i n ( 2n ) J } 


(A7 ) 


* 

Recalling  that  0  =  g—  we  obtain,  to  lowest  order  in  e. 


8*^  +  $<%>2iw1W  FK-(kY)§lr  (tjf) 


where  £  =  £0+kz  .  Integrating  this  one  obtains 


e  =  e  +  y  (t-r)  +  •  .  • 


(A8) 
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where  the  neglected  terms  are  nonsecular  and  of  order  e,  or  secular 

and  higher  order.  Finally,  noting  that  a  and  the  oscillatory  term  tend  to 

zero  as  e-*0,  (A6)  may  be  inverted  to  yield 

t  -  t  -  .  (A9) 

Afio(I-l) 

This  may  be  inserted  in  (A4),  (A7)  and  (A8)  yielding  the  orbits  needed 
to  perform  the  integration  in  (7). 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


FIGURE  CAPTIONS 


Normalized  density  D  as  a  function  of  the  distance  z  (in  cm) 
measured  from  the  anode  along  the  drift  tube.  Wiggler 
wavelength  2tt/I<  =  3  cm,  y  =  1.68,  guide  field  H0  =  2.2  kOe, 
e=6H/H0  =  1/1024,  resonance  harmonic  number  1=2,  radius 

of  anode  R0  =  1.2  cm.  Observation  point  for  lower  curve  has 

-2 

coordinates  (x,y)  =  (0,10”  )  cm;  for  upper  curve  (x,y)  = 

(0,0.5)  cm. 

Relationship  between  anode  and  typical  gyro-orbits  at  observation 
point  0  with  coordinates  (x,y)  =  (0,0.5)  cm. 

Normalized  density  D  as  a  function  of  the  distance  z  along 
the  drift  tube  when  e  =  1/128.  Other  parameters  are  the 
same  as  those  in  Fig.  1. 

Normalized  density  D  as  a  function  of  the  distance  z  along 

the  drift  tube  when  e  =  1/64.  The  observation  point  has 

_  o 

coordinates  (x,y)  =  (0,10  )  cm.  Other  parameters  are  the 

same  as  those  in  Fig.  1. 
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Fig.  5 


Fig.  6 


Normalized  parallel  temperature 
2  2 

T  =  (ck/&|n0|)  (yTz/mc  )  as  a  function  of  the  distance 
z  along  the  drift  tube.  For  the  upper  curve  e  =  1/128; 
the  lower  curve  has  e  =  1/1024.  The  observation  points 
for  both  curves  have  coordinates  (x,y)  =  (0,0.5)  cm. 
Other  parameters  are  the  same  as  those  in  Fig.  1. 

Normalized  density  D  as  a  function  of  the  distance  z 
along  the  drift  tube  when  y  =  2.17,  e  =  1/128  and  the 
resonance  harmonic  number  i  =  3.  Other  parameters  are 
the  same  as  those  in  Fig.  1. 
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